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' Abstract 

We obtain explicit upper and lower bounds on the norms of the spectral 
projections of the non-self-adjoint harmonic oscillator. Some of our results 
apply to a variety of other families of orthogonal polynomials. 

\o ' 

^ ■ 1 Introduction 

o 

m 

O I We consider polynomials p„ which are orthogonal with respect to a complex weight 

^ ■ a on [0, oo) in the following sense. We suppose that p„ is of degree n and 

(•no 

Pm{x)Pn{x)a{xf dx = 5m,n 



J3 



^ . for all non-negative integers m, n. (All of our statements and proofs can be rewrit- 

ten with (0, oo) replaced by R, and we will not keep repeating this point.) If a > 
and Pm are real-valued, then they are orthonormal in L^((0, oo), (t(x)^ dx) in the 
usual sense, but for complex-valued a such an interpretation is not possible. Our 
goal is to obtain bounds on the quantities 



oo 



Nn= \pn{x)(j{x)\ dx 

Jo 

for all n. 

This problem arose in the context of the non-self-adjoint harmonic oscillator 

{Hf){x) = -nx) + z'x'f{x) 
acting in i^^(R) for some complex z. In this situation the relevant weight is 

a{x) = e-^'"'/2 
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and Nn is the norm of the spectral projection P„ of H associated with its nth 
eigenvalue, A„ = z'^{2n+ 1). In the numerical literature Nn is called the condition 
number of the eigenvalue A„. Numerical calculations in indicated that \\Pn\\ 
increases at an exponential rate a.s n —>■ oo, and it was proved in [Ij that there was 
no polynomial bound on ||Pn|| for this and certain other Schrodinger operators. 
The super-polynomial rate of increase of the associated resolvent norms in the 
semi-classical limit was proved in |2] by a method which was greatly generalized in 
jH]. For certain classes of operators with analytic coefficients it was recently proved 
that the resolvent norms increase at an exponential rate in the semi classical limit, 
However, the precise exponential constants have not been identified in any 
example. 

A consequence of our theorems is that there exists a positive critical constant tz 
such that the 'spectral expansion' 

oo 
n=0 

is norm convergent ii t > tz and divergent if < t < t^. Our method provides 
explicit upper and lower bounds on tz but not its precise value. 

The problem may be reformulated as finding the norms of (f)n{x) = Pn{x)(T{x) 
in L^((0, oo), dx), where 0^ are obtained by applying a modified Gram-Schmidt 
orthogonalization process to the functions x"'a{x). This procedure is modified in 
the sense that we require 

/■oo 

/ (j)m{x)(l)n{x) dx = 5m,n 

Jo 

without any complex conjugates. This is equivalent to requiring that (pm and 
0* (x) = 0„(x) form a biorthogonal system in L'^{{0, oo), dx) in the sense that 

for all m, n. If P„ is the (non-orthogonal) projection 

Pnf = {fA*n)4>n 

then PmPn = 5m,nPn for all m, n and it is easily seen that 

llPnII = Nn. 

In order to make some progress with this problem, we make the following assump- 
tions on the weight a. We assume that (j{z) is an analytic function of z in the 
sector S = {z : \ arg(2;)| < a}, and that it is positive on the real axis. We also 
assume that _ 

/•oo 

/ x"|a(e*^x)|2 dx < oo 
Jo 
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for all n > and 16*1 < a, in order that p„ should be well-defined. Our most 
important condition is that 



for all 1^1 < a and all r > 0, where > and < sq < 1. Our main theorem 
provides a lower bound on iV„ for the weight x — >■ a{e^^x) under these assumptions. 
Examples of such weights are given in Section 3. Finally in Section 4 we compare 
the bounds obtained with numerical evidence. 



2 The Lower Bound 

Let {pn}'^=o denote the standard orthonormal sequence of real- valued polynomials 
with respect to the positive weight o"^ on (0, oo). We define 

where z E S and x > 0. If z > then 



by making the change of variable zx = u. By analytic continuation the same holds 
for all complex z & S. We are interested in obtaining a lower bound on the quantity 



(2) 




n,z 




for complex z E S. Note that Nn,z = 1 for all positive real z. 



Theorem 1 Under the assumption ^ we have 



Nn,z ^ CfjSg 



2n-l 



(3) 



provided z = re and \6\ < a. 



Proof We have 



n,z 




Now 



n-l 



Pn{zx/sgr) = z'^Sg'^r >„(x) + kjPj{x) 

j=0 
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for constants kj which we need not evaluate. By the orthogonahty of the polyno- 
mials, we have 



poo 



Jo 



n-1 

\pn{zx/sgr)a{x)\'^ dx = 

3=0 

-2n 



> 

The statement of the theorem follows. 
We next consider the example 

aiz) = z^l\-'' 
where 7 > -1 and /3 > 0. If r > and |^^| < 7r/(2/?) then 

|a(re'^) | = r^/^g-r^ cos(0/3) _ cea{ser) (4) 

where sq = {cos{6j3)y^'^ and cg = Sg"'^'^. After replacing (0, 00) by (—00,00), the 
particular choice 7 = and (3 = 2 leads one to the study of the Hermite polynomials 
with a complex scaling, which is relevant to the non-self-adjoint harmonic oscillator. 
The choice (3 = 1 leads to the Laguerre polynomials L^. As far as we know, all 
other choices lead to non-classical polynomials. 

The following theorem provides a more general type of weight satisfying and 
can itself easily be generalized. 

Theorem 2 // 

n 

a{x) = exp{— ^ CjX^} 
i=i 

for all X G (0, 00), where Cj G R for all j and Cn > 0, then a satisfies (0j provided 
\e\ < TT/{2n). 

Proof We have to find kg > and sg G (0, 1) such that 



Cj cos{j9)r^ < kg + Cjslr-" 



^3 

3=1 3=1 



for all r > and \6\ < 7r/(2r7,). The validity of such an inequality depends upon 
the coefficient of r". We achieve the required bound cos{n6) < < 1 by putting 



Sg = {(l + cos(n^))/2} 



l/n 



Note If (0, 00) is replaced by R in the above theorem, we must also assume that 
n is even. 
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3 The Upper Bound 



It is surprisingly difficult to obtain an upper bound on A^„, and we treat only two 
cases. We start with the orthonormal sequence of Laguerre polynomials, associated 
with the weight a{x) = e~^/^ on (0, oo). We have 

n 
r=0 



where 



satisfies 

by virtue of the general inequality 



(r!)2(n — r)! 

bn,r\ < 27^! 

+ < 2"+V!s! 



The following theorem provides an upper bound on Nn,z which complements the 
lower bound of Theorem ^ 

Theorem 3 If a{x) = e~^'/^ and z = e*^ then 

for all n > 0, provided \6\ < 7i/2 and sg = cos{6). 
Proof We start with the equality 

Nn,z = Sq'^ / \pnie'^x/se)a{x)\'^ dx, 



which is proved as in Theorem^ We have cg = 1 and sg = cos{6) by (^. We 
deduce that 



^ /"OO " ^ 

— ^9 / y \bn,rbn,s\s0 ^ G dx 

•^0 r,s=0 

~ Jo r\s\ 



r,s=0 

n 



T I el 
/ n \ ^ 

r I 



\r=0 

^ -2n-lr)4n+2 

^ Sn Z 
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Note This proof can be extended to more general weights provided suitable bounds 
on the coefficients bn,r can be obtained, but in general this is not easy. 

We next consider the non-self-adjoint harmonic oscillator. The orthonormal se- 
quence of polynomials corresponding to the weight a{x) — 1"^ is given by 
= knHn{x), where 

fc„ = 7r-i/^2-"/2(n!)-^/2 
and Hn are the Hermite polynomials 

We will need the following lemma. 

Lemma 4 If r, s are non-negative integers then 

/oo „ 
^2r+2s^-x^ dx < 7rV22'-+V! S\ 
-oo 

Proof The left hand-side equals 

r<x> 

/ u'^+'e-'s-^/^ ds = r(r + s + l/2) 
< 7r^/^r(r + s + 1) 

In the following theorem we restrict attention to the case of even integers; the 
treatment of odd integers is very similar. 

Theorem 5 Let z = e*^ where \9\ < ti/A, and put sq = (cos(2^^))^/^. Then 

A^2n,.<7r(n + l)V224"+V"-^ 
for all non-negative integers n. 

Proof We start with the identity 

n 
r=0 

where 

(-l)"-'-22'-~"^(2n)! 
" 7rV4(n-r)!(2r)! ' 
In the following chain of inequalities we will use 

2-^-^ VF+T < < 2-^V^(^ + 1) 
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for all non-negative integers r; this is proved using induction and Stirling's formula. 
Following the method of Theorem 121 we have 

/oo 



r,s=0 
n 

r,s=0 



^ (ra - r)! (2r)! 



n! fr!)2\2 



[n — r)\ r! (2r)!^ 



\r=0 



..=0 (20!, 

n 



< Ss^^''-\n + 1)-V222n j^^^ 277r(r + 1) 



4 The Spectral Expansion 

Let H denote the non-self-adjoint harmonic oscillator acting in L'^(R,) , with eigen- 
values Xn = z'^{2n + 1) and spectral projections If the right hand-side of the 
expansion 

oo 

e-^* = E e"'"*^n (5) 

n=0 

is norm convergent, then by comparing the action of the two sides on the eigen- 
functions 0„ we see that they coincide on a dense subspace, and hence on the whole 
of L\K). 

If we put 

= limsupn~Mog(||P„||) 

n— >oo 

then our theorems imply that < < oo provided < \6\ < tt/4. They also 
provide explicit upper and lower bounds on s^- 

Theorem 6 The spectral expansion (0j is norm convergent ift > tz = s^j (2 cos(20)) 
and is norm divergent if < t < t^. 
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Proof For t > tz the terms of the series decrease at an exponential rate, while for 
Q <t <tz they are not uniformly bounded in norm. 

5 Numerical Results 

The non-self-adjoint harmonic oscillator (P) has eigenvalues A„ = z'^{2n + 1) and 
eigenfunctions 

for n = 0,1, where kn are normalization constants, Hn are the Hermite polyno- 
mials, and I arg(2)| < 7r/4. 

Theorem 7 // P„ is the nth spectral projection of H and z = re*^ then 

lirninf n-Mog(||P„||) > log(sec(2^)). 

Proof This follows directly from Theorem^ upon observing that ||-Pra|| = N^^z and 
se = cos(20)V2. 

We have previously evaluated these norms numerically for z^ = c = y/i, i.e. 9 = 
Tc/ 16. See k^^-* in Table 4 of p7. It appears from the computations there that 

lim n-Hog(\\Pj) ~ 0.40 

which is considerably larger than the lower bound 0.079 of Theorem [7| 

We now report on a more systematic numerical investigation of the spectral pro- 
jections of (H)). We evaluated crn{9) = y^||P„|| / ||Pn-2|| for various n and 9 using 
Maple. (This was easier than evaluating ||Pn|| / ||-Pn-i|| because different algorithms 
are needed for even and odd n.) The method used was the same as that described 
in [H sect. 4.3]. We put Digits := 200, and included enough terms of the sequence 
determining the eigenvector to achieve stability. For each 9 it appeared that o"„(^) 
was an increasing function of n, so the limiting value is probably larger than the 
computed value. For 9 = the operator H is self-adjoint, and the projections have 
norm 1. As stated earlier one must restrict 9 to the range |^| < 7r/4. The results are 
shown for n = 100 in Table 1. The second column lists the constants ^ = sec(2^) 
(rounded down) associated with the lower bound of Theorem ^ The fourth column 
lists the constants 4s^^ = 4sec(26') (rounded up) associated with the upper bound 
of Theorem|Sl The final column lists the values of ^{9) = exp(tan(2^)), for reasons 
explained below. 

The approximations fi{9) were obtained by the following non-rigorous method. 
For even values of n the eigenfunction 0„ of H is an even function of x which is 
concentrated around the points ±xo, where Xq is defined below. On the positive 
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9 sf aioo(^) Asf M 





0.025 
0.05 
0.1 
0.15 
0.20 



1 

1.012 
1.051 
1.236 
1.701 
3.236 



1 

1.165 
1.369 
1.953 
3.062 
6.282 

Table 1 



4 

4.050 
4.206 
4.945 
6.806 
12.945 



1 

1.172 
1.384 
2.068 
3.961 
21.708 



half-line the semi-classical analysis of ^, Sect. 2] suggests that for large enough 
T] > 

is an approximate eigenvector of H with approximate eigenvalue A, where Xq = t], 
ipi = irj, %Ij2 = —iz*, and A = (l + z*^)?]^; in the notation of we are putting c = 
and a = 1, and are ignoring the term involving ip^. 

If is a positive integer and we put r] = {n/ cos(2^)}^/^, then a direct calculation 
shows that A = 2nz'^, which equals the nth eigenvalue of H to leading order as 
n — >• oo. This suggests that 

/q- + Xo)pdg 

" ^ |/o~0(. + xo)2ds| 

|/r^e-2'/'i^-'^2^' ds\ 
= exp{ntan(26')}. 

In view of the crude character of the approximations above, the similarity of crioo(^) 
and iJ,{9) in Table 1 is interesting. We conjecture that a more detailed semiclassical 
analysis might yield the correct asymptotic constant. This also seems the best hope 
for treating more general non-self-adjoint Schrodinger operators. 
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